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Abstract 

■ We extend the A;_l -factorization formalism for exclusive photoproduction of vector mesons to 
the production of electroweak bosons. Predictions for the jp Z^p and pp — > ppZ^ reactions 

CN . are given using an unintegrated gluon distribution tested against deep inelastic data. We present 

$Ih ! distributions in the Z^ rapidity, transverse momentum of Z^ as well as in relative azimuthal angle 

between outgoing protons. The contributions of different flavours are discussed. Absorption effects 

■ lower the cross section by a factor of 1.5-2, depending on the Z-boson rapidity. We also discuss 
the production of Z^ bosons in central inclusive production. Here rapidity and (xip^i, xip^2) distri- 

(~| . butions of Z^ are calculated. The corresponding cross section is about three orders of magnitude 

larger than that for the purely exclusive process. 
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I. INTRODUCTION 



There has been recently much experimental progress in the field of central exclusive 
production. The observation of exclusive central dijets py], as well as charmonia//i"'"/i~- 
pairs[2| has clearly demonstrated the feasibilty of detecting exclusive final states at collider 
energies. These are good prospects for possible future studies at the LHC, addressing a wide 
range of physics problems, from Higgs physics, to the investigations of the QCD-Pomeron 
and hadronic structure of the produced particles. For recent reviews see for example ^] for 
theory /phenomenology, and [4] for experiment. 

The mechanism of the reaction strongly depends on the centrally produced particle, in 
particular its spin, parity, C-parity and internal stucture. For heavy vector quarkonia such as 
J/\E' and T the photon-pomeron fusion is the dominant mechanism (for recent calculations 
see e.g. [s], [6|). The same is expected for the gauge boson 0, Isj. The dominant reaction 
mechanism is shown in Fig{TJ 



hi hi 




FIG. 1: Diagrammatic representation of the photon-Pomeron fusion mechanism of exclusive Z - 
boson production in hadronic collisions. 



Here, the essential ingredient is the 7p Z^p subprocess, which proceeds through the 
qq component of the virtual photon. There is a strong similarity to the production of qq 
vector mesons, and one would expect the QCD description of this process to follow from 
the color dipole/fc^-factorization approach to vector meson production (for a review, see 
(sl) by straightforward modifications [7.,^]. The important distinction to the case of vector 
meson production is the fact that the qq pair coupling to the Z^-boson can be put on its 
mass-shell. In the impact parameter space color-dipole formulation, this requires to continue 
the light-cone wavefunction of the Z° to the region of complex arguments. The resulting 
highly oscillatory integrands are however not straightforward to handle 0]. 

In the momentum space representation given in this work, the situation is more trans- 
parent, and the numerics poses no special problems. 

Previous calculations of the pp ppZ^ process made use of the equivalent-photon ap- 
proximation (EPA) and did not include absorption effects. In the EPA only total cross 
section or rapidity distribution of the Z^ boson can be calculated. In this paper we use the 
formalism of jl, H| to perform the calculation at the amplitude-level, which allows us to cal- 
culate other differential observables (e.g. in Z° boson transverse momentum or correlation 
in relative azimuthal angle between outgoing protons) and to include absorption corrections. 

The cross sections for exclusive Z^ production for both the Tevatron and LHC are very 
small. In fact a recent search for exclusive Z'^p^] only puts rather generous bounds on the 
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cross section. The exclusive events are characterized by large rapidity gaps between centrally 
produced Z° bosons and very forward or very backward final state nucleons. Another 
process with these features is the inclusive double-diffractive production of Z^, which to 
our knowledge was not previously calculated the literature ^ . In the latter case the in 
the central rapidity region is associated with low-multiplicity hadronic activity. 



II. 7p ^ Z^p PHOTOPRODUCTION PROCESS 
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FIG. 2: A pQCD diagram for the jp — > Z^p amplitude at large 7p cm-energy. The gluon exchange 
ladder is modelled by the unintegrated gluon distribution of the target. 



Before we go to the hadronic reaction let us start from the real photoproduction process 
depicted in FiglJl The forward production amplitude can be written very much the same 
way as for exclusive photoproduction of vector quarkonia (see for example j^): 

M{W,A^ = 0) = W^ [\z [d^k- ^ . (2.1) 

As was already shown in 0,0], only the vectorial part of the Z'^gg-coupling contributes to 
the forward amplitude, 

sin2evK 

is the relevant weak vector coupling, I^f is the weak isospin of a quark of fiavour /, ej is its 
charge, and Qw is the Weinberg angle. The imaginary part of Af can be obtained from the 
results given for vector meson production with the 7^ vertex in [9|. Performing azimuthal 
integrals one obtains [gJ: 



Jo ^ 

(2.3) 



^ So far only the single-difFractive contribution was estimated in the literature, see e.g. 11 1 
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with 



Aof{z, e) = m) , e) = [z' + (1 - z^T^^ , 

K ~r III J 



and 



^ ^ fc^ + / A;^ — ml — 

„2 Z\ -| Jill J 



The strong couphng enters at the scale = ma.x{k^ + m'j, n^}. The unintegrated gluon 
distribution T{x,h?) is sampled at a; = Cskewed ■ M^/W^ with Cskewed = 0.41. This shifted 
value of X simulates the prescription of [12] to obtain the skewed distribution from the 
diagonal one, and is valid for the particular gluon distribution we use j9|. Now, let 



(2.4) 



then, the integration domain z & [0, 1] must be split into [0, 2;_]U[2;_, 2;4.]U[2;+, 1]. Apparently, 
within the subdomain z^] the denominator in Eq. fl2.ip can go to zero, which means that 
the qq state after the interaction can go on-shell. This leads to a rotation of the complex 
phase of the dominantly imaginary amplitude. For z G one must use the Plemelj- 

Sokhocki formula 

^ ^^172 r7W+ ^^S{k'' + m}-z{l-z)Ml), 




k^ + m}- z{l - z)Ml -ie k^ + m) - z{l - z)Ml 

(2.5) 

where PV denotes the principal value integral. It can be evaluated as 



'0 k —T Jo k —T J(i+A)r2 k —T 

+ Af{z,T') log(A) 



(2.6) 



for an arbitrary positive value of A. Here = z{l — z)M'^ — m'j is positive in the relevant 
integration domain. 

Another distinction in comparison to the vector- meson (VM) photoproduction is worth a 
comment. Effectively, we replace the non-perturbative light cone wave-function of the VM 
by the propagator: 
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While in the case of vector-mesons, the hght-cone wave-function will suppress the endpoint- 
region 2;, 1 — z -C 1, no such suppression of asymmetric qq pairs is available here. Incidentally, 
in impact parameter space asymmetric pairs correspond to large dipole sizes [13], and it 
is precisely the wave-function suppresion of large dipoles, which leads to the dipole-size 
scanning property jl, [13] of VM production amplitudes. Therefore, strictly speaking, the 
Z° production cross section is not purely perturbatively calculable, but one must rely on 
the ability of the color-dipole//cx-factorization approaches to properly factorize the large 
dipole/ infrared contributions. Compare this to the scaling contribution of large dipoles to 



the transverse DIS structure function Ft{x,Q ) = 2xFi{x,Q ) [13|, |l5j at large Q 



Finally, we note, that we restore the real part of the amplitude by substituting 

Af = {i + p)QmAf , (2.8) 

where p = tan(7rAip/2), and Aip = dlog'^mAf/dlogW'^, and the 7p —>■ Z^p amplitude 
within the diffraction cone is given by 

M{W\ A^) = M{W\ A^ = 0) exp[Bt] , (2.9) 

where t = — A^ and the running diffraction slope is taken as 

B = Bo + 2a',ff log (w^VW^o ) , (2-10) 

with Bo = 3.5 GeV-^, a'^^j = 0.164 GeV"^ and Wo = 95 GeV 0. 



III. pp ppZ^ EXCLUSIVE HADROPRODUCTION 

Assuming only helicity conserving processes the Born amplitude for the pp ppZ^ 
reaction is a sum of amplitudes of the processes shown in Fig{T] and can be written through 
the amplitudes of photoproduction processes 7/^2 — > -Z'°/i2 or 7/12 — > Z'^h2, discussed above, 
in the form of the vector 

^h!h2-.h^h2Z<^iPi'P2) = ei V4:7iaemFi (ti) ^ ■ M^*h2^z"h2 (Q?; S2, t2) 

ZiZi 

+ e2V4:TraemFi{t2)-^ ■ M^*h,~,z"hAQl^Si,ti) . (3.1) 

Above Qi = —ti and Q2 = — ^2 are virtualities of photons, Fi is the familiar Dirac elec- 
tromagnetic form factor of the proton/antiproton and Pi, P2 are transverse momenta of 
outgoing protons. In the present analysis we shall use a simple parametrization of the Dirac 
electromagnetic form factor Fi taken from Ref . [l7| . 

The dependence of the the virtual photoproduction subprocess amplitude on space-like 
virtuality of the photon is in practice entirely negligible (recall that -C M|). 

In the hadroproduction process one has to include additional absorption corrections. The 
relevant formalism for the calculation of amplitudes and cross-sections was reviewed in some 
detail in Ref. [5|. Here we give only a the main formulas. The basic mechanisms are shown 
in FigEi 
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FIG. 3: Diagrammatic representation of the exclusive production amplitude with inclusion of 
absorptive corrections. The absorptive corrections are calculable in terms of the hadron-hadron 
elastic S-matrix. 



Inclusion of absorptive corrections (the 'elastic rescattering') leads in momentum space 
to the full, absorbed amplitude 

Mip^,p,) = J ^S,iik)M^'\p,-k,p, + k)=M^'\p„p,)-6Mip„p,). 

(3.2) 



With 



S^iik) = (2vr)25(2)(fc) - iT(fc) , T(fc) = a?!,is) exp ( - ^B^^k') , (3.3) 



where at Tevatron energy y/s = 1800 GeV crff^ (s) = 76 mb, Bel = 17 GeV-^ (l8| , the 
absorptive correction 6M reads 

6M{p,,p,)= J ^y^T{k)M^'\p,-k,p, + k). (3.4) 

The differential cross section is given in terms of M as 

= J M\^ dydtidt2d(P , (3.5) 

where y is the rapidity of the Z'^-boson, ti, t2 are four-momentum transfers squared, and 
is the angle between transverse momenta Pi and Pj- 



IV. INCLUSIVE DOUBLE DIFFACTIVE PRODUCTION OF Z' 



The purely exclusive process discussed so far can be characterized by large rapidity gaps 
between the centrally produced and forward/backward emitted protons/antiprotons. Can 
the rapidity gap method be used to identify this process? The discussed exclusive process is 
not the only one with double rapidity gaps. The inclusive double-pomeron cross section for 
Z° boson production was apparently not previously calculated. The mechanism is depicted 
in FigH 
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FIG. 4: Pomeron-Pomeron fusion representation of the central diffractive inclusive Z production. 
The hard subprocess is viewed as a collision of partons of the Pomeron. 



Following Ingelman and Schlein pJi] , one may try to estimate the hard diffractive process 
by assuming that the Pomeron has a well defined partonic structure, and that the hard 
process takes place in a Pomeron-Pomeron collision. Then the rapidity distribution of Z^- 
bosons would be obtained from 

da(pp PpZ'^X) ■r-^ _ r, \ { Df 2\ -D/ 2\ \\ 

= K (^{QfQf Z ; X1X2S) {xiqj {xi, fx ) X2qf {x2, /i ) + (1 ^ 2) 1 . 

(4.1) 

Here 

Mz y Mz _„ 

xi = —^e', X2 



Cqfqf^zo{s) is the elementary "zeroth-order" flavour- dependent cross sections (see e.g. (20|). 
K in Eg. (14. II) stands for the Drell-Yan type i^'-factor which includes approximately pQCD 
NLO corrections [20| . 

The effective 'diffractive' quark distribution of flavour / is given by a convolution of the 
flux of Pomerons /ip(xip) and the parton distribution in a Pomeron g//ip(/5, /i^): 

dxjpd/35{x - xjp(3)qf/-sp{/3,iJ,'^) fjp{xjp) = / /ip(a;ip)g//ip( — , . 

(4.2) 

The flux of Pomerons fjp {xjp ) enters in the form integrated over four-momentum transfer 

Mxjp)= r^^^ dtf{xjp,t), (4.3) 

with tmin,tmax being the kinematic boundaries. 

Both pomeron flux factors /ip(xip, t) as well as quark/antiquark distributions in pomeron 
were taken from the HI collaboration analysis of diffractive structure function and/or from 
the analysis of diffractive dijets at HERA 2l|. The factorization scale is taken as /x^^ = M|. 
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Besides the Pomeron-exchange, one must also include the secondary Reggeon-exchange 
contribution, which dominates at larger xip. Consequently, a number of interference con- 
tributions arise, which are shown diagramatically in Fig. [51 As we wish to use the results 
of the HI analysis of diffractive DIS at HERA, we have to omit a number of interference 
terms, which would involve Reggeon-Pomeron interference structure functions that have 
been neglected in the HI analysis. Incidentally, within pQCD, there is no reason why inter- 
ference terms should be small, and in fact, they enter the diffractive structure functions in 
the maximal possible way [2^. While one obtains good fits of HERA data, even omitting 
the IR — IP interference, the rather unphysical values of the so-extracted Reggeon trajectory 
parameters show the limitations of such a procedure. Clearly though, a full reanalysis of 
Hl-data is not warranted for our purpose of obtaining a first estimate of the Z°-production 
cross section. 

It is however obvious, that the naive factorization prescription of Eq. fl4.ip cannot be 
correct. It neglects rescattering effects of incoming protons shown in FigEl Indeed such 
diagrams quantify the probability that protons emerge intact out of the interactions region 



in a regime where the typical events are highly inelastic and many channels are open [23 



Here we restrict ourselves to only a 'minimal' scenario of factorization breaking induced 



by eikonalized multiple scatterings, following the formalism of Ref.[2J] (for a somewhat 
modernized version, see (25|). 

We do not enter here the debate on the possible relevance of multiple-Pomeron vertices 
(see for example the reviews [1, 2^ and references therein), but keep in mind that our 



treatment of absorption may require a revision after better knowledge on soft interactions 
at the LHC has been acquired. 

The relevant formulas are most easily written in impact parameter space. As in practice 
ti ~ — 9i , where is the transverse momentum transfer to proton i, we can write 

/ip(xip,t) = /ip(xip,q) = /ip(xip)exp[-5qr^] , (4.4) 

where B = B^xjp) = Bn + 2a^ log(l/xip) is the xip-dependent diffractive slope. We 



follow the ifl-analysis [21|, and use the central values of their fit B^) = 5.5 GeV ^,aip = 
0.06 GeV~^ Then, in impact parameter space, we have 

/ip(xip, b) = J /ip(a;iP, q) exp[-ibq] = ^^^^^^Ph^] = fjpi^Tp) ^ip(a;ip, b) . 

(4.5) 

Here 

/ip(xip) = fjp{xp,q = 0) = j <fbfjp{xjp, b) . (4.6) 
Now, we should make the following replacement in the cross section: 

/ip(a;ip,i)/ip(a;ip,2) = j d^bd'^bi£b25^ib - bi + 62) /ip(a:;ip,i, &i)/ip(xip,2, 62) 

^ J d^bd^b,d^b25\b-b, + b2)Sl,^{b)fjp{xjp,MfA^JP,2M)- 

The treatment of absorptive corrections is in fact fully analogous to the one required for 
77 collisions in heavy ion collisions, compare e.g. Eq.(2.2) in Ref. 27 . 
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(e) 

FIG. 5: Mueller-Kancheli diagrams for the process pp ^ p + gap + Z^X + gap + p. Interference 
diagrams of type (c,d), which would involve IP — IR interference structure functions are neglected. 
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FIG. 6: Absorptive correction to the IP-exchange contribution. 



TABLE I: Survival probabihty factor for purely exclusive and inclusive double diffractive processes. 



v^[GeV] 


pp ^ p + Z^X + p 


pp^Y^ + ZoX + Y2 


1960 


~ 0.05 


~ 0.06 


14 000 


~ 0.025 


~ 0.03 



Using the form fl4.5l) of the Pomeron-flux, we observe, that the Born-level cross section 
will be multiplied by the effective survival probability factor 

S^(xip,i, X]p,2) = j d'^bd^bi£b2Sl^,{b)5^{b - 6i + 62)^1? (a;iP,i, bi)tjp{xjp^2, ^2) 

= . B / r / d'bSlMexp 

47r5i2(Xip,i,Xip,2) J 

Here Bi2{xjp^i, x-sp^2) = B{x-sp^i)B{x-sp^2)/ {B{x-sp^i) + B{x-sp^2))- In fact, due to the very small 
Pomeron Regge-slope in the Hl-fit the xip^i — Xjp 2 dependence can be safely neglected. A 
two-channel model for the absorption factor S'^^^(6), is described in the appendix. It yields 
the numbers given in Table [H 

V. RESULTS 

Before we go to hadronic reactions let us first present predictions for the 7p Z^p 
reaction. In Figd we show the total cross section as a function of photon-proton center of 
mass energy W. In this calculation we have used the unintegrated gluon distribution from 
Ref. [l5| and the slope parameter B taken from Eq. fl2.10l) . The cross section grows quickly 
with the energy from the kinematical threshold Wth = niz + mp. At typical HERA energy 
W = 200 GeV the cross section is of the order of 10~^ nb, i.e. too small to be measured. 
However, it grows quickly with energy and at = 10 TeV it is already of the order of 1 pb. 
We show not only the cross section with the full amplitude (including all flavours) but also 
results with three (u+d+s: dotted), four (u+d+s+c: dashed) and five (u+d+s+c+b: solid) 



4^nf.TTD 1 . XtD r,) . 



(4.7) 
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flavours. At low energies it is enough to include only light flavours, while at large energies 
all flavours must be included. 




FIG. 7: The total cross section for — > Z p reaction as a function of photon-proton center- 
of-mass energy for the Ivanov-Nikolaev UGDF. The dotted line includes: u+d+s, the dashed line 
u+d+s+c and the solid line: u+d+s+c+b. 



In FigE] we show distributions in rapidity for the pp ppZ^ (Tevatron) and pp — *• 
ppZ^ (LHC) without (black thin solid) and with (grey thick solid) absorption effects. The 
Born approximation cross section calculated here is much larger than that calculated in 
the dipole approach in Refs. 0, [1] . Generally the absorption effects lower the cross section. 
The effect depends on the rapidity. Absorptive corrections for exclusive Z^ production 
are bigger than for the exclusive production of J/\l/ jl] and T [6|. This is due to the 
fact that for heavy particle production on average higher four momentum transfers (and 
hence less peripheral collisions) are involved than for lighter particles. Analogously as for 
photoproduction 7j9 — Z^p in Figj9] we show the distribution in Z-boson rapidity in the 
Born approximation calculated with different number of flavours included in the calculation. 
While for the Tevatron energy it is enough to include four flavours (u,d,s,c) at the LHC 
energy flve flavours must be included. At LHC the inclusion of the b quarks increases the 
cross section by about 20 %. 

In Fig{TD]we show separate contributions of photon-pomeron and pomeron-photon fusion 
mechanisms as well as the sum of both processes. We wish to stress the fact that in rapidity 
distributions all interference phenomena dissapear if absorption corrections are neglected. 
At LHC the two contributions are better separated which leads to the camel-like shape with 
minimum of the cross section at j/ ~ 0. 

The cross sections at the Tevatron energy i/i = 1.96 TeV is rather smah. Recent search 
for the exclusive Z° production [13] has found only upper limit for this process. There is 
a hope that at the LHC it could be measurable. One should remember, however, that in 
practice one can measure either fi'^fi~ or e+e~ pairs. This means one can expect a sizeable 
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(a) (b) 

FIG. 8: Rapidity distribution of the exclusively produced for the Tevatron (left) and LHC 
(right) energies. The black thin solid line corresponds to Born amplitudes and the grey thick solid 
line (red on-line) includes absorption effects. 




(a) (b) 

FIG. 9: The influence of different flavours on rapidity distributions of for Tevatron (left) and 
LHC (right) energies. The dotted line includes: u+d+s, the dashed line u+d+s+c and the solid 
line: u+d+s+c+b. 



background from the 7*7* l^l^ processes lOj ^ • In order to get rid of this type of the 
background some cuts could be helpful. 

In Fig{TT] we show transverse momentum distribution of the exclusively produced Z^. 
The distribution peaks at pt ~ 0.3 GeV and extends to relatively large transverse momenta. 
This is in clear contrast to the photon-photon processes where the corresponding transverse 
momenta of the lepton pair would peak at much lower transverse momenta. Imposing a 
lower cut on the lepton pair transverse momenta would cut off the unwanted photon-photon 
background. 

There are definite plans that both ATLAS and CMS main detectors will be supplemented 



^ Within the Standard Model, the 77 transition is absent at the tree level. In fact, single-Z" 

production in 77 collisions has up to now not been observed experimentally. 
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FIG. 10: The photon-pomeron and pomeron-photon contributions for the Tevatron (left) and LHC 
(right) energies. No absorptive corrections are included here. 




0.5 1 1.5 2 0,5 I 1.5 2 

Pi (GeV) p, (GeV) 



FIG. 11: The transverse momentum distribution of Z'^ bosons at ^fs = 1960 GeV (left) and y/s = 
14 TeV (right). 

by several forward detectors. In principle, having two forward detectors could allow to 
measure two outgoing protons in coincidence. This could allow studying correlations between 
outgoing protons. As an example in FigH^] we show the distribution in relative azimuthal 
angle between outgoing protons. Quite different distributions are obtained for the Tevatron 
{pp collisions) and LHC {pp collisions). This effect is of the interference nature and was 
already discussed for exclusive J/\l/ production [5|. In contrast to the reaction with Z*^ boson 
the relative azimuthal angle distribution for the photon-photon processes peaks sharply at 
~ 180° [13] ■ Therefore imposing extra cuts in the azimuthal angle space should further 
diminish the photon-photon background opening a chance to measure for the first time the 
exclusive Z^ production in proton-proton collisions. 

Finally let us present our estimate of the inclusive double-diffractive contribution of FigJH 
In FigHSlwe show the cross section with pomeron exchanges only (dashed) and with both 
pomeron and reggeon exchanges included (solid). This cross sections have to be multi- 
plied in addition by the gap survival probabilities from Table 1. In this calculation xj^"^ 
= 0.1 was assumed. This means cuts on longitudinal momentum fractions of outgoing pro- 
tons/antiprotons. Even after including absorption corrections the inclusive double-pomeron 
contribution is a few orders of magnitude larger than the purely exclusive cross section. The 
rapidity distributions from Fig{T3] are more narrow than the purely exclusive distributions 
shown earlier. This is partially related to the cuts on longitudinal momentum fractions and 
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f (deg) 



V (deg) 



FIG. 12: Distribution in relative azimuthal angle between outgoing protons at = 1960 GeV 
(left) and x/s = 14 TeV (right). 

is of purely kinematic origin. The cross section for inclusive double pomeron contribution 
is certainly measurable at LHC and is bigger than the cross section for single-diffractive 
production at Tevatron [28^. 
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FIG. 13: Rapidity distributions of inclusive double-pomeron production of boson for the Teva- 
tron (left) and LHC (right). The solid lines include both pomeron and reggeon exchanges while 
the dashed lines only pomeron exchanges. In this calculation the fit (a) from Ref. 21[ was used. 
No absorptive corrections were included here. 



In Fig{Tl]we show two-dimensional distribution in pomeron/reggeon momentum fractions 
(a:]P^i, a:ip^2)- The large mass of the Z^ boson causes that the small values of \ and 2 
are not accessible kinematically. This is more evident for the Tevatron energy. This is also 
the reason for much smaller cross section for double- diffract ive production of Z^ for Tevatron 
compared to LHC. 
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FIG. 14: Two dimensional distributions in (xjp^i, xip^i) for the Tevatron energy W = 1960 GeV 
(upper panels) and for the LHC energy W = 14000 GeV (lower panels). The left panels include 
only pomeron exchanges while the right panels both pomeron and reggeon exchanges. In this 
calculation the fit (a) from Ref. 21[ was used. No absorbtion corrections were included here. 



VI. CONCLUSIONS 

We have extended the fc^-factorization to the exclusive production of bosons. The 
production amplitude was calculated using an unintegrated gluon distribution [ISi] adjusted 
to inclusive deep inelastic structure functions. The so obtained 7p — > Zp amplitude served 
as an input for the evaluation of the pp ppZ^ process. Our results obtained with bare 
(i.e. without absorption) amplitudes are by a factor of 3 larger than those obtained earlier 
in the dipole approach. We have analyzed the role of individual flavours. For low energy it 
is enough to include only light flavours while at high energies all flavours must be included. 

Compared to earlier works in the literature we have taken into account absorption effects. 
The absorption effects depend on the Z-boson rapidity and lower the cross section by a factor 
of 1.5-2. As for the exclusive T production the larger rapidity, the larger the absorption 
effect. 
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Very small cross sections are obtained both for Tevatron and LHC This means that 
possible background must be studied. The is measured via e~^e~ or decay channels. 



Recently the CDF collaboration [10| presented a first estimate of the upper limit for the 
exclusive Z° production at the Tevatron. Their limit is about three orders of magnitude 
larger than our predictions. This demonstrates difficulties to measure the exclusive process. 
The situation will improve at LHC (larger cross section, larger luminosity), but even there 
it will be rather difficult to measure the cross section. 

In a detailed analysis a background from the 77 l^l^ and 7IP ^ 7* ^ 
(sub)processes must be taken into account. We have found that distributions in (lepton 
pair) transverse momentum as well as in relative azimuthal angle between outgoing protons 
can be very useful to separate out the background 77 processes. 

To our knowledge for the first time in the literature, we have estimated inclusive double 
diffractive production of Z^ using diffractive parton distributions obtained recently from 
the analysis of the proton diffractive structure functions/diffractive dijets performed by the 
HI collaboration at HERA. We have calculated cross section assuming Regge factorization 
as well as inculding absorption effects leading to the factorization breaking. Rather large 
inclusive double-diffractive cross sections were found at LHC. A future experiments with 
forward instrumentation of main LHC detectors (ATLAS and ALICE) should provide new 
results concerning hard diffraction. This will allow to further investigate the mechanism of 
Regge-factorization breaking observed already for soft total single and double diffraction. 
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VIII. APPENDIX: TWO-CHANNEL MODEL FOR THE ABSORPTIVE COR- 
RECTIONS 

Here we present the details of the two-channel model used for the evaluation of absorptive 
corrections to central inclusive Z°-production. It improves over a single-channel desription 
by taking into account some inelastic shadowing corrections. As physical states, one would 
include the proton \p) and some effective low-mass states |A^*) with proton-quantum num- 
bers (representing e.g. resonances and the diffractively excited Ntt, Ntttt, ... components. 
The multichannel eikonal will then be an operator acting on the tensor-product space of 
physical states \NiNj) = \Ni) \Nj), where Ni = p,N*. The proton-proton S'-matrix in 
this space is now given by 

S{b) = exp^-^u{s,b)g®gj, (8.1) 
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where the opacity v{s^h) is given by 



So 



(8.2) 



and we stick to an oversimphfied model, in which all matrix elements have the same h 
dependence given by Tjp{b). The 6-space profile was taken in the Gaussian form: 



Tip (6) 



1 



exp 



27r(5o + 2a' \og{s/s^)) " L2(5o + 2a' log(s/so)) 



(8.3) 



The values of the bare pomeron parameters g^^, Aip as well as the parametrisation of Tjp{b) 
are given below. Notice that below, we do not distinguish between protons and antipro- 
tons, and will always refer to proton-proton scattering even when discussing results for the 
Tevatron. 

In the two-channel case, where all inelastic excitations are subsumed in a single effective 
state |A^*), the matrix g is written as 



1 + 5 7 
7 1-6 

It has the eigenvalues 

Ai,2 = i± v'^^T^, 

and the physical states can be expanded into S-matrix eigenstates as 



•4) 



(8.5) 



(8.6) 



Now we turn to the gap survival probabilty and evaluate the effective Sl^f^^ib) which enters 
Eg. (14. 71) . We distinguish different final states: 

A. pp ^ p + Z^X + p 

First let us the case of the proton-proton final state. Here we need to substitute 



Slbsib) — ^ {pplig ® 9)S{b)\pp) = {pp\{g^g)ew[-7;^is,b)g(g)g]\pp) 



(8.7) 



To evaluate the matrix element {pp\ . . . \pp), we should expand the protons into eigenstates 
of the S'-matrix according to Eq.( ]8.6p : 



{pp\ ® g)s{b)\pp) = ICfl'Xle""'"/' + IC^^l^Xle-"'"/' + 2 \Cf\^\Cl\^X,X2e-'''' 



-Z/A1A2/2 



(8.8) 



Here we suppressed the arguments of z/ = z/(s, b). 
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B. pv^Yi + Z^X + Y2 

If protons in the final state cannot be measured, we need to sum over all excitations 
^1,2 £ {p^N*}, and we should substitute 

slsib) — $^|mi^2|(^®^)^(b)M|' = (H(^®^)'^'WM 

= {pp\{g®g) exp[-z/(s,b)^(g)^]|pp) 

= |Cf I^A^e-"^? + \Cl\^\1^e-'''-'^ + 2 \Cl\^\Cl\\\i\2fe-''^'^' . (8.9) 



Equivalent equations can be found in [25[, who we largely follow in choosing 7 = 0.55, 6 = 0. 
Then 

|p> = i=(|l) + |2)),|iV> = ^(|l)-|2>), (8.10) 

with eigenvalues A12 = (1 ±7). The bare Pomeron parameters used in the parametrisation 
of the opacity (18. 2p with Sq = 1 GeV^, are taken as 

glp = 27mb, Aip = 0.11, Bq = 9 GeV-^ a' = 0.14 GeV-^ . (8.11) 

These parameters are so adjusted, that we obtain reasonable values for the total cross section 
atot, the elastic cross section aei, as well as the elastic slope Bgi. They are obtained from 



^tot = 2 I d'bm ,a^, = l d'bV'ib) , fie. = \^^^^ 



(8.12) 



where the impact-parameter space forward amplitude r(6) is given by 

T{b) = l-{pp\S{h)\pp). (8.13) 

For the energy of Tevatron Run I, y/s = 1800 GeV, we obtain atot = 78.5 mb, = 16.7 mb, 
and Bel = 17.2 GeV'^ For the LHC energy of ^/s = 14TeV, this oversimplified model 
predicts atot = 106 mb, aei = 26 mb, and Bei = 19.8 GeV~^. 
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